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Abstract 

Drought is a prolonged reduction in precipitation which may result in significant effects on local 
ecosystems, watersheds, and water sources, with the potential to have negative impacts on human 
populations. Normalized indexes can be used to characterize drought magnitudes and frequencies 
using historical data records on temperature, rainfall, snowpack, soil moisture, streamflow, and/or 
other water supply indicators. This study examined the applicability of two drought indices: 
Standardized Precipitation Index (SPI) and Streamflow Drought Index (SDI). The SPI and SDI 
estimate the magnitude of observed rainfall (SPI) and discharge (SDI) deviations from the long¬ 
term average. 

In this study, SPI and SDI were used to determine droughts from 1970-2000 in 7 North Georgia 
catchments. Because other studies have shown the importance of time on index performance, SPI 
and SDI were calculated on data aggregated for 1, 3, 6, 9, 12, and 24 month intervals. The 
performance of the indices was evaluated using the strength of the relationship between SPI and 
observed discharge and between SPI and SDI using linear regression. The goodness of fit was 
measured by a statistic referred to as the coefficient of determination (r 2 ). The r 2 is the amount of 
variation in the dependent variable explained by the variation in the independent variable. Two- 
way Analysis of Variance (ANOVA) models were conducted to determine how watershed and 
temporal scale (independent variables) influenced the SPI-discharge and the SPI-SDI regression 
relationships. The dependent variable in these models was the coefficients of determination 
estimated using the linear regressions. 

Hydrological drought severity was analyzed on an annual basis and revealed the following trends 
for SPI and SDI: 1) The timing of severe drought events varied slightly across watersheds; 2) the 
most extreme drought events occurred in 1978, 1981, 1986-1988, and 1999-2000; and 3) 1986- 
1988 events were similar in severity to the 1999-2000 drought events. Watersheds varied 
temporally in their hydrologic response to drought, which may be a function of localized recharge 
or anthropogenic factors. Based on i^s analysis, discharge and SPI had a more robust relationship 
with longer time steps (9 - 24 months); SPI and SDI relationships revealed an analogous pattern. 
Regions deficient in discharge data can use SPI as a viable alternative to SDI; whereas, SDI can 
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function as a proxy for SPI in catchments lacking precipitation data at timescales exceeding six 
months. 

Keywords: Drought, Standardized Precipitation Index (SPI), Streamflow Drought Index (SDI), 
discharge, streamflow, Georgia 
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Introduction 

One of the most subtle, yet devastating weather catastrophes is drought. Drought is the 
reduction in rainfall relative to normal conditions (Tannehill, 1947). Because drought intensity, 
duration, and frequency are stochastic in nature, there exists little consensus on a standard 
definition for what constitutes drought conditions (Tannehill, 1947; Wilhite and Glantz, 1985; 
Quiring et al., 2007). However, scholars concur on four common indicators of drought conditions: 
meteorological (e.g., reduction in average precipitation), agricultural (e.g., reduced soil moisture 
resulting in crop damage), socio-economic (e.g., shortages of drinking water), and hydrological 
(e.g., declining reservoir or river levels) (Byun and Wilhite, 1999; Wilhite and Buchanan-Smith, 
2005; Nalbantis and Tsakiris, 2009). Drought, regardless of how it is assessed, can have profound 
effects on social, economic, and ecological systems (Dahm et al., 2003; Bernard et al., 2003). 

The United States (US), despite its well-developed infrastructure for accessing and distributing 
water, has experienced severe droughts. In 1998, drought-related wild fires in Florida burned over 
1920 square kilometers resulting in more than $500 million in damages. That same year Canada 
had the highest fire occurrence season in a quarter century, and Mexico suffered the worst drought 
in seventy years resulting disaster zone declarations (NOAA, 2013). Droughts during 1998-2002 
and 2006-2009 in the Southeastern US caused severe crop damage, disrupted hydro-electric power 
generation, and resulted in water shortages. These events resulted in seven to ten billion dollars in 
drought damages during 1998 alone (Quiring et al., 2007; Pederson et al., 2012). Repeated water 
shortages in the headwaters of the Alabama-Coosa-Tallapoosa (ACT) and Apalachicola- 
Chattahoochee-Flint (ACF) river basins have placed Georgia in the middle of a tri-state water 
conflict (Chapman and Peck, 1997). This legal conflict over scarce water resources underscores 
the importance of having effective tools for identifying drought conditions so that coordinated and 
effective water resource management can be implemented to mitigate the socio-economic impacts 
of drought (Singh et al., 2015). 

Because drought has such important socio-economic ramifications, there have been efforts to 
develop mathematical models, such as indices, for assessing drought (Wilhite et al., 2007). 
Generally these drought indices are designed to quantify intensity, duration, and/or spatial extent 
of precipitation deficits (Mishra and Singh, 2010). These indices have been used to assess drought 
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effects on lake levels, river discharge, and groundwater recharge (Vicente-Serrano et al., 2012b). 
Normalized drought indices have been proposed to measure extreme weather events because they 
are inexpensive, can apply to various time scales, are relatively simple to calculate and understand, 
and require only one type of data (e.g., precipitation, Rimkus et al., 2013). These indices measure 
the magnitude and direction of the deviation from the norm of the long-term average (precipitation, 
discharge, or other water supply indicators). The Standardized Precipitation Index (SPI) was 
originally formulated to calculate meteorologic drought conditions using rainfall records (McKee 
et al., 1993); while, the Streamflow Drought Index (SDI) is a modified version of SPI that uses 
discharge data (Nalbantis, 2008). Because the SPI can be calculated across various time scales 
(e.g., 1-24 months), it has been used to assess drought effects on reservoir storage, stream flow, 
and irrigation (Giddings et al., 2005; Quiring et al., 2007). 

These two indices differ in the types of droughts they assess most effectively. SPI is 
particularly valuable for identifying meteorological drought conditions, but does not factor in 
geologic characteristics of the study region (McKee et al., 1995). Mountainous topography can 
lead to seasonal drought variability through fluxes of rainfall and snowmelt (Vasiliades et al., 
2011) and can, under certain circumstances, influence SPI-derived drought estimates. Conversely, 
SDI integrates more effectively the geologic characteristics that can confound studies of 
hydrological drought that use SPI, because SDI accounts for water table fluctuations that are 
controlled by rainfall, surface runoff, and subsurface flows. SDI has been broadly applied as 
illustrated by its use to delineate hydrological drought in the Athens, Greece Metropolitan 
Watershed (Nalbantis, 2008), Spanish Ebro River Basin (Vicente-Serrano et al., 2012), 
northwestern Iran (Tabari et al., 2013), and the European Neman River Basin (Rimkus et al., 2013). 
Given their practicality for use in real-world applications, it is surprising that so few studies (e.x. 
Nalbantis and Tsakiris, 2009 and Al-Faraj et al., 2014) directly compare the effectiveness of these 
two indices to delineate hydrological drought events. A comparison of SPI and SDI could be used 
to quantify regional historical hydrological drought formation, duration, and intensity or assess if 
these indices can act as proxies for basins lacking sufficient meteorological or hydrological data. 

The purpose of this study was to evaluate the efficacy of the SPI and SDI for recording 
hydrologic drought. There were two specific goals of this research: 1) to determine how watershed 
(N=7) and temporal scale (1, 3, 6, 9, 12, and 24 month(s)) influence the relationship between 
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streamflows and SPI and 2) to assess how watershed and temporal scale affect the relationship 
between SPI and SDI. A comparison of SPI and SDI across seven North Georgia watersheds will 
provide a better understanding of the usefulness of both indices for quantifying hydrological 
drought. 

Methods 

Study Areas 

The study area covers a 60,000 km 2 region that includes parts of Georgia, North Carolina, 
and South Carolina in the headwaters of the Alabama-Coosa-Tallapoosa (ACT) and Apalachicola- 
Chattahoochee-Flint (ACF) river basins (Fig. 1). These basins are critically important, because 
they are the source of drinking water for millions of people and the subject of ongoing disputes 
about water use allocations (Chapman and Peck, 1997). This region has abundant rainfall, with 
an annual average of 1000 to 2000 mm, typical of humid subtropical climatic regions (Konrad, 
1997; Diem and Mote, 2005; Gotvald et al., 2009). The region is comprised of the Piedmont and 
Blue Ridge geologic provinces and is characterized by mountainous terrain underlain by igneous 
and metamorphic aquifers (Gotvald et al., 2009). To enhance the generality of this study, seven 
watersheds (Table 1) were selected because they showed similar geology (e.g. crystalline bedrock), 
watershed size (400-1000 km 2 ), and available hydrologic data (1970-2000) from the National 
Oceanic and Atmospheric Administration (NOAA) and US Geological Survey (USGS). Appendix 
A contains more detailed descriptions of each watershed. 

Given limitations regarding the availability of long-term precipitation and discharge data, 
this study focused on watersheds having at least thirty years of discharge (USGS) and precipitation 
(NOAA - Global Historical Climatology Network - Daily or GHCND) data. Thirty years of data 
are recommended for SPI to calculate an accurate estimate of average climate conditions (McKee 
et al., 1993); however, 60 or more years would be necessary to account for the cyclical nature of 
droughts (Wu et al., 2005). Discharge and precipitation records from 1970 to 2000 were analyzed 
for this study because these data were readily available and records were relatively complete. The 
stations used as the source of the precipitation data were selected based on the following 
conditions: 1) located either within the stream’s watershed or 30 km buffer of the USGS station, 
2) no data gap longer than three sequential months, and 3) a total of no more than 18 months of 
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missing data. Since SPI and SDI were assessed on monthly timescales, all daily rainfall and 
discharge were aggregated into cumulative monthly totals before analysis. 
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Figure 1. Study region depicting state boundaries, regional topography, seven study watersheds, 
USGS stream gauges (stars), and NOAA precipitation stations (open squares). 

























Table 1. Watershed name, site number, size, weather station name, and missing data for the 
seven watersheds used to calculate the SPI. 
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Watershed Name 

Site 

Number 

Size 

(Sq Km) 

Weather Station 

Missing Data 

Chattahoochee River 

1 

822 

Cornelia 

None 

Chestatee River 

2 

393 

Cleveland 

None 

Sweetwater Creek 

3 

616 

Dallas 

None 

Coosawattee River 

4 

613 

Ellijay 

Weather: Nov 74, Dec 75, 
Aug 79, & Feb 92 

Chattooga River 

5 

537 

Coweeta 

Weather: Oct 00 

Flint River 

6 

675 

Atlanta Hartsfield 

None 

Middle Oconee River 

7 

1021 

Gainesville 

None 


SPI & SDI Analysis 

SPI calculates the difference between precipitation over a defined time period (Pi) and the 
period of record mean precipitation (P m ) scaled by the variation over the period of record (i.e. - 
standard deviation [SDp]) (McKee et al., 1995; Heim, 2002). 


SDp 


SPI is a dimensionless number where negative values indicate drought and positive values indicate 
wet periods (Bordi et al., 2001; Bonaccorso et al., 2003). This model assumes the underlying 
rainfall data are normally distributed (Panu and Sharma, 2002). Efforts should be made to use 
these models on complete data sets, because incomplete data sets may not represent a normal 
distribution (Mishra and Singh, 2010). 

One of the goals for this study was to analyze USGS discharge records using a modified 
SPI analysis, referred here as the Streamflow Drought Index (SDI, developed by Nalbantis and 
Tsakiris, 2009). More recently, SDI has also been called the “Standardized Streamflow Index” or 
SSI (Vicente-Serrano et al., 2012a; 2012b; Tsakiris, et al., 2013). SDI could serve as an alternate 
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drought indicator because discharge data are readily available and contain more complete records 
than many of NOAA precipitation data. SDI was calculated as 

_ Qi-Qm 

SDI = - 

sd q 

where Qi is the discharge over a defined time period, Q m is the period of record mean 
discharge, and SDq is the discharge’s standard deviation over the period of record. 

SPI and SDI data were calculated using SPI gamma distribution software 
(http://drought.unl.edu/MonitoringTools/DownloadableSPIProgram.aspx) available through the 
National Drought Mitigation Center (NDMC). To investigate how time scale influences model 
behavior, both SPI and SDI drought indices were calculated on 1, 3, 6, 9, 12 and 24 month 
timescales. Missing data were converted to zero as suggested by the NDMC programming 
instructions (2016) since less than two percent of data was absent. A modified Mann-Kendall Test 
(Mann, 1945; Kendall, 1955; Hamed and Rao, 1998) was used to show that missing values do not 
affect long-term weather data analysis if less than seven percent of data is involved (e.x. Sen Roy 
and Balling, 2004; Zhang et al., 2011). 

GIS Analysis 

Site selection required three objectives to be met: 1) find a USGS stream gauge with 
discharge data, 2) develop a watershed for the USGS stream gauge, and 3) locate respective 
weather stations. First, an Excel file containing decimal degree coordinate locations of USGS 
streamflow gauges was converted to point shapefiles using the “Add XY Data” function. X was 
defined as longitude and Y as latitude. This same process was used to develop the weather station 
point layer. Once the streamflow point layers were imported into ArcMap, a 30 km buffer was 
created around the streamflow point layer using the “Buffer” Geoprocessing Tool. Study 
watersheds were created by merging a HUC12 basin layer and then reshaping the polygon to end 
at the streamflow gauge through visual analysis of a one meter DEM and surface flow lines. 
NOAA weather stations were then selected either within the 30 km buffer or study watershed 
boundary. 



7 


A series of descriptive maps were created to understand how land use, geology, ecoregions, 
and dam locations may have influenced each watershed. A 2001 raster land use file (USGS) was 
clipped to the individual study watershed features, where it was converted into a vector file. The 
attribute tables of each watershed’s land use was then exported into an Excel (Microsoft 
Corporation ® 2010) spreadsheet and the count of land use values was converted to percentages 
by dividing the count by the total number of pixel counts. Geology and ecoregion polygon layers 
were imported into ArcMap from external sources. These were then clipped to the individual 
watershed layers. Once clipped, a field was added to the attribute table to calculate area (in square 
meters) for each feature within the layer using the calculate geometry tool. The attribute tables of 
each watershed’s geology and ecoregions were then exported into an Excel spreadsheet and 
percentages were calculated by dividing the individual layer areas by the total study area. Dam 
layers were created by importing the Georgia National Inventory of Dams (NID) point layer into 
ArcMap and clipping it the study areas. Since the Chattooga River (Site #5) also spanned North 
and South Carolina (NC and SC), analysis of an aerial image was conducted to determine other 
dam locations. These were then verified by the Corp of Engineers of NC and SC NID database. 

ArcMap 10.1 (ERSI) was used to conduct all GIS analysis. No ArcMap extensions were 
used. Excel files, shapefiles, and raster layers were downloaded from the Internet (Table 11, 
Appendix D), and imported into ArcMap. Land use, geology, and dam location maps for each 
watershed are included in Appendix A. 

Statistical Analysis 

The performance of the indices was evaluated using the strength of the relationship 
between SPI and observed discharge and between SPI and SDI. The strength of these relationships 
was assessed using coefficients of determination derived from linear regression models (Sokal and 
Rohlf, 1995). Least-squares linear regression is a statistical technique used to quantify the slope 
and intercept for the linear relationship between independent (x) and dependent (y) variables. The 
goodness of fit is measured by a statistic referred to as the coefficient of determination (r 2 ). The r 2 
is the amount of variation in the dependent variable explained by the variation in the independent 
variable (Sokal and Rohlf, 1995; Mittlbock and Heinzl, 2002; Morid et al., 2006). An r 2 of 1 
indicates 100 % of the variance in y is explained by x. For this study, least-squared linear 
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regressions relating SPI (independent variable) to discharge (dependent variable) and SPI 
(independent variable) and SDI (dependent variable) were performed at each time scale (N = 6) 
and for each watershed (N = 7) combination (i.e., N = 42 total) using Microsoft Excel (ver. 2013). 
Due to the large number of regression models computed, the threshold probability for statistical 
significance was lowered to p = 0.001 (i.e., 0.05/42) using a Bonferroni correction (Sokal and 
Rohlf, 1995; Armstrong, 2014). 

Two-way Analysis of Variance (ANOVA) models were used to determine how watershed 
and temporal scale (independent variables) influenced the SPI-discharge and the SPI-SDI 
regression relationships. The dependent variable in these models was the coefficient of 
determination estimated using the linear regressions. Two-way ANOVAs were selected for this 
analysis, because it could be used to assess both watershed and temporal scales simultaneously. 
The interaction term (site * time scale) was omitted due to a lack of replication. Tukey HSD was 
used because it corrects for the number of pairwise comparisons (Sokal and Rohlf, 1995). Both 
Two-way ANOVA statistical models and comparisons were performed in IBM SPSS Statistics 
software (ver. 21, IBM Corp., 2012). 

Results and Discussion 

Characteristics of the Study Area 

Draining an average of 668 ± 203 km 2 (mean ± SD), the seven study watersheds are located in 
two Level III ecoregions (Griffith et al., 2001). These watersheds are located primarily within the 
Piedmont ecoregion with approximately one third of the area falling within the Blue Ridge 
ecoregion (Table 2). Only two rivers, Chattooga and Coosawattee, are located completely within 
the Blue Ridge ecoregion; while Sweetwater Creek, Flint River, and Middle Oconee Rivers occupy 
the Piedmont ecoregion. Ecoregions are defined by their biotic (e.g., vegetation and wildlife) and 
abiotic (e.g., geology, physiography, climate, soils, and land use) characteristics (Omemik, 1987; 
1995). These characteristics can affect a watershed’s hydrology (Kokkoken et al., 2003). 
Dominated by gneiss, schist and granite, the Piedmont ecoregion is considered hillier than the 
mountainous Blue Ridge ecoregion (Griffith et al., 2001). The Blue Ridge ecoregion is a rugged 
terrain mix of igneous (mafic and ultramafic), metamorphic (low and high grade), and sedimentary 
(bouldery colluvium) geology (Griffith et al., 2001). Land uses (2001) in the study areas were 
predominately forested with less than 30% developed and agriculturally cultivated (Table 2). 
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Other land uses in the watersheds included wetlands and barren landscapes, but these uses 
composed less than 5% of the land area (Table 2). All watersheds had dendritic drainage systems 
normally associated with the hilly metamorphic and igneous terrain of the Appalachian Mountains 
(Chapman and Peck, 1997). 


Table 2. Mean ( x ) and standard deviation (SD) of percentages of watershed area by level III 
ecoregion and land use category (2001) for the seven study areas. 

Watershed Site Number 



1 

2 

3 

4 

5 

6 

7 

X 

SD 

Ecoregion 

Piedmont [45] 

72.4 

69.6 

100.0 

0.0 

0.0 

100.0 

100.0 

63.1 

45.0 

Blue Ridge [66] 

27.6 

30.4 

0.0 

100.0 

100.0 

0.0 

0.0 

36.9 

45.0 

Land Use 

Barren 

0.3 

0.5 

0.5 

0.1 

0.1 

0.5 

1.1 

0.4 

0.3 

Cultivated 

12.7 

9.5 

11.4 

5.1 

1.9 

13.4 

25.0 

11.3 

7.4 

Developed 

8.6 

7.3 

33.9 

6.6 

3.7 

34.9 

16.5 

15.9 

13.2 

Forested 

77.8 

82.3 

49.2 

87.9 

94.0 

42.4 

53.5 

69.6 

20.7 

Wet 

0.6 

0.4 

5.0 

0.3 

0.2 

8.8 

3.8 

2.7 

3.3 


Table 3 shows annual average basin rainfall varied 8 - 30% across the seven study 
catchments from 1970-2000. The Blue Ridge ecoregion typically experiences annual average 
precipitation over 2000 mm (Griffith et al., 2001). Located completely within the Blue Ridge 
ecoregion, the Coosawattee River (Site #4) and Chattooga River (Site #5) averaged similar annual 
precipitation. Site #5 also had the highest discharge of all of the watersheds, which may be 
connected with its land use and protection under the Wild and Scenic Rivers Act (Kent and Bayne, 
2010; Smith and Moore, 2011). Ninety-four percent of the Chattooga River Watershed is 
comprised of forested land which is protected by three national forest districts: the Highlands 
District of the Nantahala National Forest (NC), the Andrew Pickens District of the Sumter National 
Forest (SC), and the Tallulah District of the Chattahoochee National Forest (GA) (Kent and Bayne, 
2010 ). 
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Table 3. Mean annual precipitation (Pi 2) and standard deviation (Psd) among study watersheds 
including the three highest ranking 12 month SPI values and their associated drought year (3). 


Most Intense SPI Ranked Droughts (12 Month) 


Site 

Pn (mm) 

Psd (mm) 

1 st (year) 

2 nd (year) 

3 rd (year) 

1 

1519 

246 

-2.33 

(1981) 

-1.97 

(2000) 

-1.69 

(1988) 

2 

1634 

272 

-1.96 

(1978) 

-1.85 

(2000) 

-1.39 

(1988) 

3 

1368 

224 

-2.17 

(1986) 

-1.80 

(1988) 

-1.67 

(1978) 

4 

1546 

262 

-2.28 

(1987) 

-1.93 

(1993) 

-1.50 

(1988) 

5 

1819 

334 

-1.98 

(2000) 

-1.81 

(1988) 

-1.76 

(1986) 

6 

1269 

186 

-2.18 

(2000) 

-1.62 

(1999) 

-1.36 

(1986) 

7 

1383 

205 

-2.12 

(1988) 

-2.05 

(1978) 

-1.84 

(1981) 


The Flint River (Site #6) catchment was the driest basin conceivably due to the complex 
precipitation-recharge relationship (Kokkenen et ah, 2003), historically low rainfall rates (Singh, 
2015), water withdraws (Ruhl, 2005; Pulido-Velazquez et ah, 2008), or large scale dams 
(Chapman and Peck, 1997). Urban rainfall modification caused precipitation rates to vary in three 
contiguous watersheds studied within the Charlotte, North Carolina metropolitan area (Wright et 
ah, 2013). Thirty-five percent of the Flint River basin is developed which may have led to storm 
formation downwind of the urban heat islands and lower rainfall rates within the watershed. Other 
studies (Kokkenen et ah, 2003; Skoien et ah, 2003; Zelelew and Alfredsen, 2014) have suggested 
that the location of rain and/or discharge gauges may affect hydrologic modeling results; Flint 
River’s weather station is located in top of the watershed while the discharge gauge is on the 
bottom. Furthermore, the Flint’s linear watershed shape could have affected accumulation totals. 
In North Carolina, the elongated shape of Little Sugar Creek watershed may have led to lower 
accumulation of rainfall due to basin shape and orientation in relation storm motion and formation 
to three broader watersheds (Wright et al., 2013). Groundwater recharge primarily maintains 
baseflows in the southern portions of Flint River within the Coastal Plain due to low rainfall rates 
especially in the summer months (Singh, 2015). Cropland irrigation could be furthermore 
decreasing the base flow in the Flint River (Torsk et ah, 1996; Ruhl, 2005; Johnson, 2013), 
although this problem would be concentrated in the southern most portion of the Flint River 
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watershed. Dams may have also contributed decreased streamflow by increasing evaporative 
losses; Site #6 has 52 water impoundments (Table 5 and Appendix A [dam location maps]). 


Table 4. Mean annual stream discharge (Qn) and standard deviation (Qsd) among watersheds 
including the three highest ranking 12 month SDIs and their associated drought year in 
parentheses. 


Drought Severity 12 Month SDI 


Site 

Qn (m 3 /s) 

Qsd (m 3 /s) 

I s * Highest 

2 nd Highest 

3 rd Highest 


1 

8341 

2108 

-2.08 


(1988) 

-2.03 

(1986) 

-2.03 

(1981) 


2 

3835 

1064 

-2.10 


(1986) 

-1.91 

(2000) 

-1.83 

(1988) 


3 

3735 

1049 

-2.42 


(1986) 

-2.41 

(1988) 

-1.62 

(1999) 


4 

5459 

1496 

-2.50 


(1988) 

-1.96 

(1986) 

-1.82 

(2000) 


5 

7015 

1894 

-2.13 


(1981) 

-2.04 

(1988) 

-1.94 

(2000) 


6 

3588 

1105 

-2.34 


(2000) 

-2.15 

(1999) 

-1.81 

(1986) 


7 

5469 

1715 

-2.09 


(1988) 

-2.03 

(2000) 

-1.78 

(1981) 


Table 5. 

Total number of dams and their respective uses 

among watersheds. 








Watershed Site Number 




Dam Type 

1 

2 

3 

4 

5 

6 

7 

X 

SD 

Farm Pond 

0 

0 

1 

0 

0 

2 

1 

1 

l 

Fishing 

1 

0 

0 

0 

0 

0 

0 

0 

0 

Flood Control 

13 

0 

1 

15 

0 

0 

13 

6 

7 

Hydroelectric 

2 

0 

0 

0 

0 

0 

1 

0 

1 

Irrigation 

1 

0 

0 

0 

0 

0 

0 

0 

0 

Recreation 

30 

15 

45 

23 

0 

48 

18 

26 

17 

Water Supply 

0 

0 

0 

1 

5 

2 

2 

1 

2 

Other 


0 

0 

0 

1 

4 

2 

0 

1 

2 

Total 


47 

15 

47 

40 

9 

54 

35 

35 

17 
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Drought Extremes as Identified by SPI & SDI Analysis 

SPI (Table 3) and SDI (Table 4) indicated extreme drought events occurred historically in all 
seven watersheds. According to the SPI, the most extreme calendar year meteorological and 
hydrological droughts occurred during the last 22 years ofthe study (1978,1981,1986-1988,1999- 
2000). Similarly, the SDI’s highest drought severities occurred in last 19 years of the period of 
study (1981,1986,1999-2000). These results contrasted slightly from a 1943-1995 ACF and ACT 
River basins study (Chapman and Peck, 1997), which reported intense droughts from 1984-1988 
based on USGS daily mean streamflow hydrographs. Barber and Stamey (2000) analyzed Georgia 
data from 1930-2000 and found surface water shortages similar in magnitude to the 1986 drought 
from 1998-2000. These authors also found that the 1980-82 water years resulted in the lowest 
stream levels since 1954. According to my analysis, this period experienced less severe drought 
than 1986, since the twelve month SDI and SPI results revealed more severe droughts across 
Georgia during 1986 in comparison to 1980-82 drought events (Tables 3 & 4). However, at a few 
sites the SPI (Chattahoochee and Middle Oconee Rivers) and SDI (Chattahoochee, Chattooga, and 
Middle Oconee Rivers) were consistent with findings of Barber and Stamey (2000). 

SPI results varied only slightly across watersheds (Fig. 2a). The overall 12 month hydrological 
pattern indicated that droughts showed similar timing across the period of record with an intense 
drought at the end of the record. During these drought events, Georgia may be vulnerable to 
municipal, industrial, and recreational water shortages and other socio-economic impacts (Tsakiris 
et al., 2013). Bordi et al (2001) saw a decrease in wetness starting in the latter half of the 1970’s 
for Italy during a long term SPI study (1950-2000); local weather patterns or climate change may 
have been the cause (Bordi et al, 2001; Fisher et al., 2010). Another long term SPI study in China 
(Zhang et al., 2010) reported localized weather patterns and seasonal conditions caused SPI to vary 
across 42 Pearl River Basin catchments from 1960-2005. Less frequent or heavy rainfall events 
may have skewed SPI trends (Zhang et al., 2009; Fisher et al., 2010). 

Normalized indices such as SPI and SDI assume normal probability distributions and thus, 
could be sensitive to data gaps. For example, missing data, commonly treated as zeroes, register 
in models as dry events where no precipitation occurred (Quiring et al., 2007; Mishra and Singh, 
2010). From 1970-2000, all hydrologic streamflow records were complete for all study 
watersheds, so this problem was avoided in these SDI calculations. Two study sites were missing 
meteorological data, the Coosawattee River (4 months) and Chattooga River (1 month). The 
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Chattooga River (Site #5) was missing less than 1% of its weather station data (1/360 months = 
0.3%). The Coosawattee River (Site #4) had about a 1% data deficiency (4/360 months = 1.1%), 
but none of the missing data from these months occurred within the same 12 month period. 


— Cornelia —Cleveland —Dallas —Ellijay —Cowceta —Atlanta Airport —Gainesville A 




Figure 2. Twelve month drought estimates for seven watersheds (a) SPI and (b) SDI. The dashed 
line indicates the threshold (0) between wet and dry conditions. 


Other SPI studies have been challenged by data gaps. Lana et al. (2001) conducted an SPI 
study in Spain using 99 rain gauges (1961 to 1990). In this study rainfall gauges missing 2 or more 
months in a 12 month period were excluded. Sen Roy and Balling (2004) conducted a SPI study 
on 129 stations in India where 10.2% of data records were missing for all stations, but a variety of 
analyses indicated that the missing data had little effect on annual rainfall trends. Zhang et al. 
(2012) found similar results with monthly precipitation trends when they conducted a modified 
Mann-Kendall test and compared the use of the gap-fill method to replace missing precipitation 
values (1%) with the insertion of zeroes for 12 rain gauging stations. The gap-fill method consists 
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of replacing missing data values with statistically similar data from other gauging stations such as 
a nearby weather station (Zhang, 2011). Relative to these studies, the current study has relatively 
few data gaps which are unlikely to alter the study’s conclusions. 

SDI scores varied moderately across watersheds (Fig. 2b), possibly due to localized river flow 
regimes (Tabari et al., 2013). SDI-estimated extreme droughts matched recorded water shortages 
for Georgia (Quiring et al., 2007). The overall 12 month hydrological record showed a drought 
near the end of the study period. Recent studies in Iran and Iraq indicated similar drying patterns 
from 1990-2009 (Tabari et al., 2013; Rimkus et al., 2013) which may have been due to climate 
change (Al-Faraj et al., 2014) or increased water demand from metropolitan areas or agriculture 
(Barber et al., 2000). 


Table 6. Drought classification and corresponding event probabilities for 12 Month SPI and SDI 
for the seven watersheds studied. 


Drought Index value 

Category 

SPI Probability (%) 

SDI Probability (%) 

2.00 or more 

Extremely wet 

0.4 

0.0 

1.50 to 1.99 

Severely wet 

2.1 

1.3 

1.00 to 1.49 

Moderately wet 

11.1 

4.7 

-0.99 to 0.99 

Near normal 

66.1 

74.4 

-1.49 to-1.00 

Moderately dry 

11.0 

4.9 

-1.99 to-1.50 

Severely dry 

6.3 

7.7 

-2 or less 

Extremely dry 

2.6 

5.7 


From 1970 to 2000, SPI and SDI displayed similar hydrological drought severities (Table 6) 
and overall timing (Fig. 2) in North Georgia. Quiring et al. (2007) found similar patterns in Texas. 
Fourteen Pinios River Basin watersheds in Greece (Vasiliades et al., 2011) found SPI identified 
historic hydrological droughts, but were highly variable across the watersheds both in severity and 
timing due to seasonality. Rimkus et al. (2013) reported similar results in the Neman River basins 
(Europe); seasonality caused fall droughts and spring flooding which were detected using both SPI 
and SDI. In the present study, drought events identified using the SPI tended to occur prior to SDI 
droughts (Fig. 2). The Blue River Basin in Oklahoma (Liu et al., 2012) found SPI’s meteorological 
components responded more quickly in identifying severe droughts than a hydrological drought 
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index, Standardized Runoff Index. The Reconnaissance Drought Index (a meteorological drought 
index) proactively predicted water availability trends in the Diyala watershed spanning the Iraq 
and Iran border prior to SDI (a hydrological drought index) (Al-Faraj et al., 2014). Hydrological 
drought components tend to respond slower than meteorological components, which may be due 
to the complexity of the precipitation-runoff relationship (Quiring et al., 2007; Liu et al., 2012; 
Vicente-Serrano et al., 2014). Streamflow may be a better indicator for hydrological drought since 
it takes into account the multiple factors affecting regional river systems such as: water allocations 
and withdrawals (Pulido-Velazquez et al., 2008; Johnson, 2013), stream network structure (Skoien 
et al., 2003; Lopez-Moreno et al., 2013), infiltration rates, and topography (Hewlett and Hibbert, 
1967; Kokkenen et al,. 2003). 


Factors Influencing the Relationship between SPI and Stream Discharge 

To assess the SPI’s reliability as a drought indicator, the strength of the linear association (as 
measured by coefficient of determination or r 2 ) between SPI and measured discharge were 
analyzed (Fig. 3). SPI and discharge were linearly related with r 2 values averaging 0.48 ± 0.19 (x 
± SD) across all basins and timescales. SPI versus streamflow r 2 values deviated from the average 
about 1 - 11 % across watersheds and 2 - 90 % among temporal scales (Table 7). The use of the 
^s in further statistical analyses was supported by the fact that 97% of all regressions were 
statistically significant (Table 7). The only exception was the 24 month time interval for the Flint 
River (Site #6, r 2 = 0.345, p = 0.021) which did not meet the Bonferroni adjusted alpha of p < 
0 . 001 . 

Study site location had an important effect on the relationship between SPI and streamflow 
(Two-way ANOVA, F = 3.268, p = 0.014). Temporal factors (i.e. 1 to 24-month time scale) 
significantly influenced the SPI and discharge relationship (Two-way ANOVA, F = 56.91 1>P< 
0.0005). According to a pairwise comparison, all watersheds were similar except for Site #6 (Flint) 
and #5 (Chattooga) (Fig. 4a, Tukey HSD, p > 0.026 for all). The i^s values for 1, 3, and 6 month 
time scales were statistical similar (Fig. 4b, Tukey HSD, p > 0.178 for all), but r 2 values were 
lower than those measured at the comparable 9, 12, and 24 month time scales (Fig. 4b, Tukey 
HSD, p > 0.998 for all). The time period with the strongest linear associations between SPI and 
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discharge differed by catchment: Site 2, 3, and 6 at 9 months. Site 5 at 12 months, and Site 1, 4, 
and 7 at 24 months (Table 7). 
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Figure 3. Least-squares linear regression of SPI versus discharge rates (m 3 /s) for the Flint River 
(Site #6) on a 12 month time scale (1970 to 2000). 


Due to the time lag in the precipitation-surface runoff relationship, other hydrologic drought 
studies found SPI’s timing varied in response to discharge between study watersheds (Quiring et 
al., 2007; Vasiliades et al., 2011; Wright et al., 2013). Temporal lags between precipitation and 
streamflow can be affected by anthropogenic actions such as surface and ground water withdrawals 
for agricultural, municipal, and industrial purposes (Giddings et al., 2005; Vicente-Serrano et al., 
2014; Singh, 2015). From 1961 to 2007, 192 weather stations were analyzed with SPI to 
understand historical dry/wet patterns in Zhujiang River Basin (South China). These patterns were 
influenced by the spatial and temporal distribution of monsoons and periods of no rainfall along 
with urban land use and heat island effects (Fischer et al., 2010). A Northwest Pacific study 
(Abatzoglou et al., 2014) reported that the relationship between streamflow and SPI varied across 
21 unregulated drainage basins both seasonally and spatially, with the most robust relationship 
between SPI and discharge occurred at 6 - 10 month timescales. Six Texas watersheds showed 
similar results; SPI best correlated with discharge and SDI at 6 and 9 month timescales (Quiring 
et al., 2007). 
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Figure 4. Boxplots showing coefficients of determination (r 2 ) for the linear regressions relating (a) 
SPI and stream discharge among study sites (N = 7); (b) SPI and stream discharge at various time 
scales (N = 6); (c) SPI and SDI among study sites (N = 7); and (d) SPI and SDI at various time 
scales (N = 6). Medians are represented by the dark horizontal line. Box ends represent 25 th and 
75 th percentiles. Whiskers represent the 10 th and 90 th percentile. Gray circles represent outliers [> 
1.5 * interquartile range (IQR)]; stars represent extreme outliers [> 3 * IQR)]. Bars not sharing 
the same letters represent significant statistical difference (p < 0.001) as indicated by Tukey HSD. 
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Table 7. Measure of fit of least-squares linear regressions (r 2 ) relating SPI to observed discharge 
(Q) and SDI on 1, 3, 6, 9, 12, and 24-month time scales (1970 to 2000). Bolded values mark r 2 
values which did not met the Bonferroni corrected significance threshold of p < 0.001. 



Site 

1 Month 

3 Month 

6 Month 

9 Month 

12 Month 

24 Month 

1 

SPI v Q 

0.27 

0.33 

0.36 

0.63 

0.75 

0.82 


SPI v SDI 

0.39 

0.51 

0.61 

0.69 

0.77 

0.82 

2 

SPI v Q 

0.26 

0.30 

0.34 

0.67 

0.61 

0.58 


SPI v SDI 

0.41 

0.52 

0.63 

0.74 

0.61 

0.59 

3 

SPI v Q 

0.28 

0.29 

0.31 

0.72 

0.70 

0.63 


SPI v SDI 

0.53 

0.56 

0.67 

0.74 

0.70 

0.65 

4 

SPI v Q 

0.20 

0.25 

0.25 

0.59 

0.61 

0.67 


SPI v SDI 

0.35 

0.45 

0.53 

0.66 

0.58 

0.67 

5 

SPI v Q 

0.26 

0.35 

0.41 

0.70 

0.80 

0.73 


SPI v SDI 

0.37 

0.53 

0.63 

0.78 

0.81 

0.76 

6 

SPI v Q 

0.26 

0.34 

0.38 

0.58 

0.55 

0.34 


SPI v SDI 

0.45 

0.60 

0.63 

0.67 

0.56 

0.34 

7 

SPI v Q 

0.26 

0.32 

0.38 

0.69 

0.67 

0.73 


SPI v SDI 

0.40 

0.53 

0.63 

0.71 

0.69 

0.71 


Factors Influencing the Relationship between SPI and SDI 

To evaluate if SDI could be used as a substitute for SPI, it was critical to establish what factors 
influence their association. In this study, the measure of fit (l^s) between SPI and SDI was 
influenced by location and timescale. SPI and SDI were linearly related with r 2 values averaging 
0.60 ±0.13 (x ± SD) across all basins and timescales. SPI and SDI r 2 values deviated from the 
average 0.4 - 11.3 % across watersheds and 0.7 - 16.3 % among temporal scales analyzed. All 
SPI and SDI regressions were statistically significant (Table 7) except the 24 month time interval 
for the Flint River (Site #6, r 2 > 0.344, p > 0.021). 

Statistical analysis of spatial and temporal factors influencing the relationship between SPI and 
SDI revealed similar results to those reported for SPI and streamflow. Watershed location had no 
effect on the relationship between SPI and SDI (Fig. 4c, Two-way ANOVA, F = 2.022, p = 0.094). 
However, time scale significantly influenced the SPI and SDI (Fig. 4d, Two-way ANOVA, F = 
14.02, p < 0.0005). The influence of time scale on the SPI and SDI relationship deviated slightly 
from its effect on the SPI and discharge patterns reported earlier. A one month interval had the 
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lowest fit for all timescales. One and three month intervals were statistically similar (Tukey HSD, 
p = 0.086) and generally lower than other periods. The 3 month time scale had intermediate fit 
values and was statistically indistinguishable from 1, 6, and 24 month time scales (Tukey HSD, p 
> 0.074 for all). Six month or greater time scales had the strongest SPI and SDI relationships and 
were statistically indistinguishable from one another (Tukey HSD, p > 0.249 for all). 

Consistent with this study, Nalbantis and Tsakiris (2009) found high correlations between SPI 
and SDI for predicting extreme hydrological events at 3, 6, 9 and 12 month time scales in the 
Evinos River basin (Greece). Recent studies found robust relationships between SPI and other 
hydrological drought indices at 9 and 12 month time scales for 8 Aragon River watersheds 
(Vicente-Serrano and Lopez-Moreno, 2005), a 24 month time scale for 22 southern Italy basins 
(Mendicino and Versace, 2006), and 3,6, and 9 month time scales in the Diyala watershed between 
Iraq and Iran (Al-Faraj et al., 2014). In contrast, the Neman catchment spanning five countries 
(Lithuania, Belorussia, Russia, Poland, and Latvia) found SPI and SDI were poorly correlated, 
most likely due to the mountainous terrain creating winter snowfalls and high spring discharge 
(Rimkus et al., 2013). 

Summary and Conclusions 

In this study, historical hydrologic droughts from 1970-2000 were evaluated using SPI and 
SDI in seven North Georgia watersheds. Hydrological drought severity was analyzed on an annual 
basis and revealed the following trends for SPI and SDI: 1) all seven catchments experienced more 
dry events than wet events based on percent probability (Table 6); 2) SDI indicated eight percent 
more normal events than extreme events based on percent probability than SPI (Table 6); 3) the 
timing of severe drought events varied only slightly across watersheds; 4) the most extreme 
drought events occurred in 1978, 1981, 1986-1988, or 1999-2000; and 5) 1986-1988 events were 
similar in severity to the 1999-2000 drought events. 

Watersheds varied temporally in their hydrologic response to drought, which may be a function 
of localized recharge or anthropogenic factors. Based this analysis, discharge and SPI had a more 
robust relationship with longer time steps (9 - 24 months), while SPI and SDI revealed an 
analogous pattern. However, it should be noted that r 2 exceeded 0.5 at the 6 month timescale for 
SDI and SDI. According to Tukey HSD analysis, catch basins were statistically similar. Future 
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studies should be conducted to determine if SPI and SDI relationships vary temporally and 
spatially in other southern ecoregions such as in the coastal plain. 

From a management standpoint, quantifying historical hydrologic drought can advance our 
understanding of drought dynamics and improve our response to these natural events. In Georgia, 
data availability is a key challenge limiting our capacity to conduct these analyses. Data gaps in 
existing gage records and an absence of gages in both USGS (discharge) and NOAA (precipitation) 
limited this study’s site selection. Future removal of stream and/or rainfall gages could limit our 
knowledge of drought, especially when 30 or more years is recommended for climate analysis. 
Accurate long-term and wide-range meteorological and streamflow monitoring networks are 
needed for water resource management especially for improved understanding of hydrological 
response patterns. Regions in Georgia’s Blue Ridge and Piedmont ecoregions that are deficient in 
discharge data can use SPI as a viable alternative to SDI. In addition, this study’s results indicate 
SDI can function as a proxy for SPI in catchments lacking precipitation data, particularly at longer 
time scales. 
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Appendix A Watershed Descriptions 

All but one of the watershed study areas are located solely in Georgia: Chattahoochee River 
Near Cornelia, Chestatee River Near Dahlonega, Sweetwater Creek Near Austell, Coosawattee 
River Near Ellijay, Chattooga River Near Clayton, Flint River Near Griffin and Middle Oconee 
River Near Athens. The Chattooga River watershed near Clayton spans Georgia, South Carolina, 
and North Carolina. Georgia has naturally abundant water resources, however, the growth of the 
Atlanta metropolitan area has sparked many debates about future water availability of Georgia. 

Georgia is a very dynamic geological state with four major geologic regions: the Valley 
and Ridge, Blue Ridge, Piedmont, and Coastal Plain. Most of the study sites were found within 
the Blue Ridge region. In northwest Georgia, the Valley and Ridge is comprised mainly of 
Paleozoic sedimentary rocks which were faulted and folded into northeast-southwest-trending 
valleys and ridges. In the northeast comer of Georgia, the Blue Ridge is comprised primarily of 
low-to-high-grade metamorphic rocks forming the North Georgia Mountains. The Piedmont is 
comprised of igneous and metamorphic rocks with lower elevations and less relief than the Blue 
Ridge. It is considered hilly and located in the middle of Georgia. The Fall Line marks the 
boundary between the Coastal Plain and Piedmont. Sedimentary rocks dating from Cretaceous to 
Cenozoic age are found within the Coastal Plain, which is located in the southern part of Georgia. 

Ecoregions denote areas with comparable ecosystems characterized by abiotic and biotic 
factors such as geology, topography, vegetation, land use, and wildlife (Omemik, 1987; 1995; 
Griffith et al., 2001). Georgia encompasses six Level III ecoregions. The Southwestern 
Appalachians ecoregion is located in the northwestern comer of Georgia, where it abuts the Ridge 
and Valley. Faulted and folded shale, sandstone, limestone, dolomite, marble, and mudstone 
dominate these two ecoregions. Rugged high-gradient outcrops of Precambrian gneiss, schist, 
quartzite, slate, conglomerate, phyllite, metagraywacke, and other metamorphosed rocks comprise 
the Blue Ridge ecoregion in the northeastern corner of Georgia. The Piedmont is a transitional 
ecoregion between the northern mountainous forest dominated ecoregions and the low gradient 
mosaic of pastures, farmland, forests, and wetlands of the southern coastal plains. 
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Figure 5. Level III Ecoregions in Georgia. 
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Figure 6. Level IV Ecoregions with watershed. 
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Table 8. The watershed land use percentages as of 2001 were primarily composed of forest and farming landscapes. The remaining 
topography was a mix of developed space, riparian environments, and barren land. 

Land Use Type Watershed Site Number 



1 

2 

3 

4 

5 

6 

7 

X 

G 

Barren Land 

0.29 

0.46 

0.46 

0.07 

0.15 

0.54 

1.12 

0.44 

0.35 

Cultivated Crops 

0.01 

0.01 

0.00 

0.17 

0.03 

0.04 

0.05 

0.04 

0.06 

Deciduous Forest 

58.76 

66.50 

23.90 

73.01 

64.25 

18.11 

38.95 

49.07 

21.97 

Developed, High Intensity 

0.10 

0.10 

1.10 

0.03 

0.00 

4.24 

0.44 

0.86 

1.54 

Developed, Low Intensity 

1.08 

0.60 

13.49 

0.91 

0.11 

12.98 

6.66 

5.12 

5.97 

Developed, Medium Intensity 

0.35 

0.23 

2.28 

0.19 

0.03 

4.74 

1.12 

1.28 

1.72 

Developed, Open Space 

7.05 

6.41 

17.05 

5.44 

3.55 

12.97 

8.32 

8.68 

4.71 

Emergent Herbaceous Wetlands 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.01 

0.00 

0.00 

Evergreen Forest 

9.89 

7.59 

20.83 

7.37 

20.87 

19.58 

5.90 

13.15 

6.92 

Hay/Pasture 

12.67 

9.50 

11.38 

4.95 

1.89 

13.31 

24.99 

11.24 

7.37 

Herbaceous 

4.51 

3.72 

3.01 

0.86 

0.36 

3.14 

7.83 

3.35 

2.48 

Mixed Forest 

3.99 

3.75 

0.45 

5.96 

7.95 

0.87 

0.59 

3.37 

2.91 

Open Water 

0.47 

0.35 

0.93 

0.12 

0.09 

1.97 

0.57 

0.64 

0.65 

Shrub/Scrub 

0.69 

0.69 

1.05 

0.72 

0.59 

0.69 

0.25 

0.67 

0.23 

Woody Wetlands 

0.15 

0.09 

4.07 

0.20 

0.14 

6.81 

3.19 

2.09 

2.66 
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Table 9. Geologic zones vary by watershed with a mix of igneous, metamorphic, and sedimentary rocks. 

Geologic Unit Watershed Site Number 



1 

2 

3 

4 

5 

6 

7 

X 

G 

Basal Lower Cambrian Clastic Rocks 

3.7 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.53 

1.40 

Catacalastic rocks 

4.4 

0.0 

0.0 

0.0 

1.7 

0.0 

0.0 

0.87 

1.69 

Felsic Orthogneiss (Granite/Gneiss) 

0.0 

0.0 

31.3 

0.0 

0.0 

0.0 

0.0 

4.47 

11.82 

Felsic Paragneiss/Schist 

0.9 

0.0 

0.0 

0.0 

0.0 

16.5 

50.1 

9.64 

18.86 

Lower Paleozoic 

3.4 

0.0 

0.0 

13.1 

0.0 

0.0 

21.1 

5.36 

8.41 

Lower Paleozoic granitic rocks 

8.7 

24.9 

31.7 

0.0 

0.0 

0.0 

0.0 

9.33 

13.48 

Mafic Paragneiss (Homblendite/Amphibolite) 

10.1 

16.3 

26.3 

0.0 

0.0 

0.0 

0.0 

7.54 

10.52 

Middle Paleozoic granitic rocks 

0.0 

0.0 

0.0 

0.0 

0.0 

4.1 

4.6 

1.24 

2.13 

Migmatite 

0.0 

0.0 

0.0 

0.0 

0.0 

79.4 

24.2 

14.80 

29.88 

Orthogneiss 

0.0 

0.0 

0.0 

0.0 

14.8 

0.0 

0.0 

2.12 

5.60 

Paragneiss and schist 

35.6 

10.6 

9.0 

0.0 

62.6 

0.0 

0.0 

16.81 

23.78 

Ultramafic rocks 

0.0 

0.0 

0.0 

0.0 

1.5 

0.0 

0.0 

0.22 

0.58 

Z sedimentary rocks 

33.2 

48.2 

1.7 

86.9 

19.4 

0.0 

0.0 

27.08 

32.24 
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Table 10. Mountainous ecoregions (Piedmont [45] and Blue Ridge [66] dominate the watersheds. 

Ecoregion Level IV Watershed Site Number 



1 

2 

3 

4 

5 

6 

7 

X 

G 

8.3.4 45a Southern Inner Piedmont 

72.4 

69.6 

100.0 

0.0 

0.0 

0.0 

29.0 

38.71 

41.70 

8.3.4 45b Southern Outer Piedmont 

0.0 

0.0 

0.0 

0.0 

0.0 

100.0 

71.0 

24.43 

42.56 

8.4.4 66d Southern Crystalline Ridges & Mountains 

27.6 

30.4 

0.0 

0.0 

100.0 

0.0 

0.0 

22.57 

36.79 

8.4.4 66g Southern Metasedimentary Mountains 

0.0 

0.0 

0.0 

99.7 

0.0 

0.0 

0.0 

14.25 

37.69 

8.4.4 66j Broad Basins 

0.0 

0.0 

0.0 

0.3 

0.0 

0.0 

0.0 

0.04 

0.10 
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Chattahoochee River Watershed near Cornelia, Georgia (Watershed Site #1) 

This river is a prominent water resource of Georgia and a key player in the Apalachicola- 
Chattahoochee-Flint River (ACF) Basin water wars. The Chattahoochee River is the largest of all 
the rivers inside the ACF Basin. Concerns over local fishing and salt water intrusion from the Gulf 
have placed restrictions on how much water can be withheld from Apalachicola Bay. Atlanta also 
uses this river as its primary water source. Approximately 430 miles in length, the Chattahoochee 
River stretches from the Blue Ridge Mountains near Brasstown Bald South to become a tributary 
of the Apalachicola River immediately north of the Georgia-Florida state boundaries, draining a 
total area of 8,770 square miles. 

Headwaters of this significant waterway were included in this study. This drainage area is 
822 square kilometers (318 square miles) and includes portions of Habersham, White, Union, and 
Towns Counties. The location of the USGS gauging station (#02331600) at latitude 34°32'26.6"N, 

longitude 83°37'21.99"W (Fig. 7) is 
referenced to North American Datum of 
1983, and lies within Habersham County, 
GA, Hydrologic Unit 03130001, on the 
downstream side of the bridge on Duncan 
Bridge Road (GA 384), 1.0 mile 

downstream from Soque River, 6.0 miles 
northwest of Cornelia, and at mile 401.4. 
The USGS discharge data period of record 
is from August 1957 to present. 

Cornelia is the respective NOAA 
weather station. It is located at 
34°30'58.783" N, longitude 83°32T.428" 
W. Its period of record spans from May 
1919 to present. 

Figure 7. USGS Gauge #02331600 for the 



Chattahoochee River (USGS, 2000). 







34 


Legend 


N 



Landuse 

1 Open Water 
Developed. Open Space 
Developed, Low Intensity 
| Developed. Medium Intensity 
| Developed. High Intensity 
| Barren Land 
Deciduous Forest 
j Evergreen Forest 
j Mixed Forest 
| Shrub/Scrub 
HI Herbaceous 
§§ Hay/Pasture 
| Cultivated Crops 
■I Woody Wetlands 

□ Chattahoochee Watershed 

□ Chestatee Watershed 
- Streams 

★ USGS Discharge Station 


0 



15 20 

Miles 


1:230,000 


Legend |\J 

Geology 



Basal Lower Cambrian Clastic Rocks 
Catacalastic rocks 
Felsic Paragneiss/Schist 
Lower Paleozoic 
Lower Paleozoic granitic rocks 
Mafic Paragneiss (Homblendite/Amphibolite) 
Paragneiss and schist 
Z sedimentary rocks 


□ Chattahoochee Watershed 

□ Chestatee Watershed 

★ USGS Discharge Station 



1:230,000 


Figure 8. Chestatee and Chattahoochee watershed land use during 2001 (top) and geology 
(bottom). 
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Chestatee River Watershed near Dahlonega, Georgia (Watershed Site #2) 

The Chestatee River flows 51 miles to become a major tributary of the Chattahoochee 
River joining it at Lake Lanier. It flows through the Blue Ridge geologic region and is part of the 
Dahlonega gold belt, which was the main target region for gold exploration during the 1829 
Georgia Gold Rush. The city of Dahlonega uses the Chestatee River for tubing, whitewater rafting, 
kayaking, and fishing. 

The headwaters of this waterway were included in this study. The drainage area is 393 
square kilometers (152 square miles) and includes portions of Lumpkin and White Counties. The 
location of the USGS gauging station (#02333500) at latitude 34°3T41" N, longitude 83°56'23" 
W (Fig. 11) are referenced to North American Datum of 1983, and lies within Lumpkin County, 
GA, Hydrologic Unit 03130001, on the left bank 250 feet upstream from Bearden Bridge on GA 
52, 2.0 miles downstream from Ballplay Creek, 2.5 miles east of Dahlonega, and 3.5 miles 
upstream from Yahoola Creek. The USGS discharge data period of record is from July 1929 to 

Cleveland (GHCND: 
USC00092006) is the respective 
NOAA weather station. It is 
located at 34° 35' 9.24" N and 
longitude 84° 45’ 56.8794" W. 
Its period of record spans from 
April 1943 to present. 

Figure 10. USGS Gauge 
#02333500 for the Chestatee 
River (USGS, 2000). 

Sweetwater Creek Watershed near Austell, Georgia (Watershed Site #3) 

This stream is the largest stream in Cobb County, Georgia. Sweetwater Creek State Park 
is named after this stream and is the site of an abandoned mill. New Manchester Manufacturing 
Company used the local waterway to construct fabric products, but in 1864, General William 
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Sherman had the mill burned to ground during the American Civil War. Its ruins stand today as a 
tourist attraction. Sweetwater Creek is a 303(d) listed stream for fecal coliform. 



Figure 11. Local mill ruins near Sweetwater Creek. Photo taken by author. 


Headwaters of this waterway were included in this study. The Drainage area is 616 square 
kilometers (238 square miles) and includes portions of Carroll, Paulding, Cobb, and Douglas 
Counties. The location of the USGS gauging station (#02337000) at latitude 33°46'35.4", longitude 
84°36'56.2" (Fig. 14) is referenced to North American Datum of 1927, and lies within Douglas 
County, GA, Hydrologic Unit 03130002, near the right bank 100.0 feet upstream from the bridge 
on Interstate 20, 400.0 feet upstream from Blair Bridge, 3.0 miles southeast of Austell, and 5.5 
miles upstream from its confluence with the Chattahoochee River. The USGS discharge data 
period of record is from May 1904 to present. 

Dallas 7 NE is the respective NOAA weather station. It is located at 33°58'59.391" N, 
longitude 84°45'7.35" W. Its period of record spans from February 1947 to present. 
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Figure 12. Sweetwater Creek flowing through Sweetwater Creek State Park. Photo taken by 
author. 

During 2001, the Sweetwater River watershed was moderately developed. Land uses were 
primarily deciduous forest, evergreen forest, developed open space, developed low intensity, and 

pasture lands. Developed land only totaled 33.93 %. 
Agricultural land totaled 11.38 %. Undeveloped 
land totaled 53.76 %. 

The Sweetwater Creek watershed is located 
in the Blue Ridge geologic region of Georgia. 
Lower Paleozoic granitic rocks, felsic orthogneiss, 
and mafic paragneiss are the main geological units. 
Paragneiss, schist, and preCambrian 
metasedimentary rocks compose the remaining 
geology of the Sweetwater watershed. 


02337000 - Sweetwater Creek near Austell, GA - March 12,1973 Figure 13. USGS Gauge #02337000 for Sweetwater 

Creek (USGS 2000). 
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Figure 14. Sweetwater watershed land use during 2001 (top) and geology (bottom) 
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Figure 15. Dam locations and totals within the Sweetwater watershed. 
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Coosawattee River Watershed near Ellijay, Georgia (Watershed Site #4) 

The Coosawattee River used to be utilized for whitewater rafting prior to the construction 
of Carter’s Lake. It starts north of Ellijay at the confluence of the Ellijay and Cartecay Rivers. 
The Coosawattee River flows through the Blue Ridge geologic region. Fishing, canoeing, and 
kayaking are common on this river. James Dickey's novel Deliverance was partially inspired from 
canoeing this river along with Chattooga and Chattahoochee Rivers. 

The headwaters of this waterway were included in this study. The drainage area is 613 
square kilometers (237 square miles) and includes portions of Fannin, Gilmer, Pickens, and 
Dawson Counties. The location of the USGS gauging station (#02380500) at latitude 34°40'30" 
N, longitude 84°30'31" W (Fig. 17) is referenced to North American Datum of 1927, and is located 
within Gilmer County, GA, Hydrologic Unit 03150102, on the right bank 0.5 miles downstream 
from State Highway 5, 2.0 miles southwest of Ellijay, and 2.2 miles downstream from the 
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confluence of the Cartecay and Ellijay Rivers. The USGS discharge data period of record is from 
October 1938 to present. 

Ellijay (GHCND: USC00093115) is the respective NOAA weather station. It is located at 
34° 40' 59.9874" N and longitude 84° 28' 59.9874" W. Its period of record spans from June 1937 
to January 2014. 

During 2001, the Coosawattee River watershed was largely undeveloped. Land uses were 
primarily deciduous forests. Developed land only totaled 6.56 %. Agricultural land totaled 5.12 
%. Undeveloped land totaled 88.48 %. 

The Coosawattee watershed is located in the Blue Ridge geologic region of Georgia. 
PreCambrian metasedimentary rocks are the main geological units. Lower Paleozoic rocks 
compose the remaining geology of the Coosawattee watershed. 



Figure 16. USGS Gauge #02380500 location in stream channel (left) and general location (right) 
for the Coosawattee River (USGS, 2000). 
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Figure 18. Dam locations and totals within the Coosawattee watershed. 


Chattooga River Watershed near Clayton, Georgia (Watershed Site #5) 

This river is recognized as a national wild and scenic river. It flows through the Blue Ridge 
geologic region and three states, starting in North Carolina, then flowing into South Carolina, and 
ending in Georgia. The Chattooga River flows 57 miles to become the main tributary of the 
Tugaloo River in Georgia. Commercially, the Chattooga is used for rafting and was featured in 
the book and film Deliverance. Also, it intersects the Ellicott Rock Wilderness. 

Headwaters of this waterway were included in this study. The drainage area is 536 square 
kilometers (207 square miles) and includes portions of Rabun (Georgia), Jackson (North Carolina), 
and Oconee (South Carolina) Counties. The location of the USGS gauging station (#02398000) 
at latitude 34°48'50" N, longitude 83°18'22" W (Fig. 20) is referenced to North American Datum 
of 1927, and is located within Oconee County, SC, Hydrologic Unit 03060102, on the left bank 
150 feet downstream from the bridge on US 76, 2.8 miles upstream from Stekoa Creek, 9.0 miles 
downstream from its confluence with Warwoman Creek, 9.0 miles upstream from its confluence 
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with Tallulah River, and 7.0 miles southeast of Clayton, Georgia. Its USGS discharge data period 
of record is from October 1939 to present. 

Coweeta Experiment Station (GHCND: USC00312102) is the respective NOAA weather 
station. It is located at 35° 4' 9.984" N and longitude 83° 26' 49.992" W. Its period of record spans 
from December 1942 to present. 



02398000 Chattooga River at Summerville, GA, March 14,1973 


During 2001, the Chattooga River 
watershed was largely undeveloped. Land 
uses were primarily deciduous forest and 
evergreen forest. Developed land only 
totaled 3.69 %. Agricultural land totaled 
1.91 %. Undeveloped land totaled 94.40 
%. 

The Chattooga watershed is located 
in the Blue Ridge geologic region of North 
Carolina, South Carolina, and Georgia. 
Paragneiss, schist, and metasedimentary 
rocks are the main geological units. 
Catacalastic, ultramafic, and orthogneiss 
rocks compose the remaining geology of 
the Chattooga watershed. 

Figure 19. USGS Gauge #0239800 for the 
Chattooga River (USGS, 2000). 
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Figure 20. Chattooga watershed land use during 2001 (top) and geology (bottom) 
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Figure 21. Dam locations and totals within the Chattooga watershed. 
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Flint River Watershed near Griffin, Georgia (Watershed Site #6) 

The FlintRiver is the second largest waterway in the ACF River Basin. The Flint River 
provides a drinking water source for much of western Georgia. Flowing approximately 344 miles 
in length, the Flint River stretches from the upper Piedmont (South of Atlanta) to become a 
tributary of the Chattahoochee River in the Gulf Coastal Plain near the southwestern Georgia state 
boundary. 

The headwaters of this significant waterway were included in this study. Its drainage area 
is 703 square kilometers (272 square miles) and includes portions of Pike, Spalding, Fayette, 
Clayton, Henry, and Fulton Counties. The location of the USGS gauging station (#02344500) at 
latitude 33°14'39", longitude 84°25'45" (Fig. 23) is referenced to North American Datum of 1983, 
and lies within Spalding County, GA, Hydrologic Unit 03130005, at the downstream side of a 
bridge pier on GA 16, 1.5 miles downstream from Shoal Creek, 5.5 miles upstream from Line 
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Creek, 10.0 miles west of Griffin, and at river mile 304.4. USGS discharge data period of record 
is from March 1937 to present. 

Atlanta Hartsfield International Airport is the respective NOAA weather station. It is 
located at 33°38'0" N, longitude 84°27'0" W. Its period of record spans from January 1930 to 
present. 

During 2001, the Flint River watershed was fairly developed. Land uses were primarily 
deciduous forest, evergreen forest, pasture lands, developed low intensity, and developed open 
space. Developed land only totaled 34.93 %. Agricultural land totaled 13.35 %. Undeveloped 
land totaled 51.71 %. 

The upper Flint River watershed is located in the Piedmont geologic region of Georgia. 
Migmatite is the main geological unit. Middle Paleozoic granitic and felsic paragneiss/schist 
compose the remaining geology of the Flint watershed. 



Figure 22. USGS Gauge #02344500 for the Flint River (USGS, 2000). 
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Figure 23. Flint watershed land use during 2001 (top) and geology (bottom) 
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Figure 24. Dam locations and totals within the Flint watershed. 

Middle Oconee River Watershed near Athens, Georgia (Watershed Site #7) 

The Middle Oconee River starts in Hall County, Georgia and flows through the Piedmont 
geologic region. The Middle Oconee River is a drinking water source for Athens, Georgia. 
Whitewater rafting and fishing are common on this river. Fecal coliform, nitrogen, and sediment 
are listed on the 303(d) Georgia Impairment List as pollutants of concern for the Middle Oconee 


River. 


The headwaters of this waterway were included in this study. Its drainage area is 1021 
square kilometers (394 square miles) and includes portions of Hall, Gwinnett, Barrow, Jackson, 
Clarke, and Oconee Counties. The location of the USGS gauging station (#02217500) at latitude 
33°56'48" N, longitude 83°25'22" W (Fig. 26) is referenced to North American Datum of 1983, 
and lies within Clarke County, GA, Hydrologic Unit 03070101, on the left bank 0.5 miles upstream 
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from US 78 and US 29 Business, 2.0 miles west of Athens, and 5.0 miles upstream from Barber 
Creek. The USGS discharge data period of record is from October 1901 to present. 

Gainesville (GHCND: USC00093621) is the respective NOAA weather station. It is 
located at 34° 17' 59.9994" N and longitude 83° 50' 59.9994" W. Its period of record spans from 
October 1891 to present. 

During 2001, the Middle Oconee River watershed was moderately developed. Land uses 
were primarily deciduous forest and pasture lands. Developed land only totaled 16.54 %. 
Agricultural land totaled 25.05 %. Undeveloped land totaled 58.41 %. 

The Middle Oconee watershed is located in the Piedmont geologic region of Georgia. 
Lower Paleozoic paragneiss, schist, migmatiterocks are the main geological units. Granitic and 
catacalastic rocks compose the remaining geology of the Middle Oconee watershed. 



02217500 - Middle Oconee River near Athens, GA - March 12,1973 

Figure 25. USGS Gauge #02217500 for the Middle Oconee River (USGS, 2000). 
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Figure 26. Middle Oconee watershed land use during 2001 (top) and geology (bottom) 
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Figure 27. Dam locations and totals within the Middle Oconee watershed. 
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Appendix B Least-Squares Linear Regression SPI vs Q Figures 


2500 -| 
2000 - 


y = 190.82X + 694.33 
r 2 = 0.2694 

p<0.0000 


A 


♦ 


1000 - 

500 - 

0 - 

-4 



~i- 1 -r 

-2 0 2 

SPI Values 


n 

4 



SPI Values 


Figure 28. Least-squares linear regression of the SPI in relation to discharge rates on (A) 1 & (B) 3 month 
time scale for the Chattahoochee River Watershed (1970 to 2000). 
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Figure 29. Least-squares linear regression of the SPI in relation to discharge rates on (A) 6 & (B) 
9 month time scale for the Chattahoochee River Watershed (1970 to 2000). 
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Figure 30. Least-squares linear regression of the SPI in relation to discharge rates on (A) 12 & (B) 
24 month time scale for the Chattahoochee River Watershed (1970 to 2000). 
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Figure 31. Least-squares linear regression of the SPI in relation to discharge rates on (A) 1 & (B) 
3 month time scale for the Chestatee River Watershed (1970 to 2000). 
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Figure 32. Least-squares linear regression of the SPI in relation to discharge rates on (A) 6 & (B) 
9 month time scale for the Chestatee River Watershed (1970 to 2000). 
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Figure 33. Least-squares linear regression of the SPI in relation to discharge rates on (A) 12 & (B) 
24 month time scale for the Chestatee River Watershed (1970 to 2000). 
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Figure 34. Least-squares linear regression of the SPI in relation to discharge rates on (A) 1 & (B) 
3 month time scale for the Sweetwater Creek Watershed (1970 to 2000). 
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Figure 35. Least-squares linear regression of the SPI in relation to discharge rates on (A) 6 & (B) 
9 month time scale for the Sweetwater Creek Watershed (1970 to 2000). 
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Figure 36. Least-squares linear regression of the SPI in relation to discharge rates on (A) 12 & 
(B) 24 month time scale for the Sweetwater Creek Watershed (1970 to 2000). 
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Figure 37. Least-squares linear regression of the SPI in relation to discharge rates on (A) 1 & (B) 
3 month time scale for the Coosawattee River Watershed (1970 to 2000). 
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Figure 38. Least-squares linear regression of the SPI in relation to discharge rates on (A) 6 & (B) 
9 month time scale for the Coosawattee River Watershed (1970 to 2000). 
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Figure 39. Least-squares linear regression of the SPI in relation to discharge rates on (A) 12 & (B) 
24 month time scale for the Coosawattee River Watershed (1970 to 2000). 
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Figure 40. Least-squares linear regression of the SPI in relation to discharge rates on (A) 1 & (B) 
3 month time scale for the Chattooga River Watershed (1970 to 2000). 
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Figure 41. Least-squares linear regression of the SPI in relation to discharge rates on (A) 6 & (B) 
9 month time scale for the Chattooga River Watershed (1970 to 2000). 
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Figure 42. Least-squares linear regression of the SPI in relation to discharge rates on (A) 12 & (B) 
24 month time scale for the Chattooga River Watershed (1970 to 2000). 
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Figure 43. Least-squares linear regression of the SPI in relation to discharge rates on (A) 1 & (B) 
3 month time scale for the Flint River Watershed (1970 to 2000). 
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Figure 44. Least-squares linear regression of the SPI in relation to discharge rates on (A) 6 & (B) 
9 month time scale for the Flint River Watershed (1970 to 2000). 
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Figure 45. Least-squares linear regression of the SPI in relation to discharge rates on (A) 12 & (B) 
24 month time scale for the Flint River Watershed (1970 to 2000). 
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Figure 46. Least-squares linear regression of the SPI in relation to discharge rates on (A) 1 & (B) 
3 month time scale for the Middle Oconee River Watershed (1970 to 2000). 
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Figure 47. Least-squares linear regression of the monthly SPI in relation to discharge rates on (A) 
6 & (B) 9 month time scale for the Middle Oconee River Watershed (1970 to 2000). 
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Figure 48. Least-squares linear regression of the SPI in relation to discharge rates on (A) 12 & (B) 
24 month time scale for the Middle Oconee River Watershed (1970 to 2000). 
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Figure 49. Least-squares linear regression of the SPI in relation to the SDI on (A) 1 & (B) 3 month 
time scale for the Chattahoochee River Watershed (1970 to 2000). 
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Figure 50. Least-squares linear regression of the SPI in relation to the SDI on (A) 6 & (B) 9 month 
time scale for the Chattahoochee River Watershed (1970 to 2000). 
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Figure 51. Least-squares linear regression of the SPI in relation to the SDI on (A) 12 (B) 24 month 
time scale for the Chattahoochee River Watershed (1970 to 2000). 
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Figure 52. Least-squares linear regression of the SPI in relation to the SDI on (A) 1 & (B) 3 month 
time scale for the Chestatee River Watershed (1970 to 2000). 
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Figure 53. Least-squares linear regression of the SPI in relation to the SDI on (A) 6 & (B) 9 month 
time scale for the Chestatee River Watershed (1970 to 2000). 
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Figure 54. Least-squares linear regression of the SPI in relation to the SDI on (A) 12 & (B) 24 
month time scale for the Chestatee River Watershed (1970 to 2000). 
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Figure 55. Least-squares linear regression of the SPI in relation to the SDI on (A) 1 & (B) 3 month 
time scale for the Sweetwater Creek Watershed (1970 to 2000). 
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Figure 56. Least-squares linear regression of the SPI in relation to the SDI on (A) 6 & (B) 9 month 
time scale for the Sweetwater Creek Watershed (1970 to 2000). 
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Figure 57. Least-squares linear regression of the SPI in relation to the SDI on (A) 12 & (B) 24 
month time scale for the Sweetwater Creek Watershed (1970 to 2000). 
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Figure 58. Least-squares linear regression of the SPI in relation to the SDI on (A) 1 & (B) 3 month 
time scale for the Coosawattee River Watershed (1970 to 2000). 
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Figure 59. Least-squares linear regression of the SPI in relation to the SDI on (A) 6 & (B) 9 month 
time scale for the Coosawattee River Watershed (1970 to 2000). 
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Figure 60. Least-squares linear regression of the SPI in relation to the SDI on (A) 12 & (B) 
month time scale for the Coosawattee River Watershed (1970 to 2000). 
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Figure 61. Least-squares linear regression of the SPI in relation to the SDI on (A) 1 & (B) 3 month 
time scale for the Chattooga River Watershed (1970 to 2000). 
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Figure 62. Least-squares linear regression of the SPI in relation to the SDI on (A) 6 & (B) 9 month 
time scale for the Chattooga River Watershed (1970 to 2000). 
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Figure 63. Least-squares linear regression of the SPI in relation to the SDI on (A) 12 & (B) 
month time scale for the Chattooga River Watershed (1970 to 2000). 
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Figure 64. Least-squares linear regression of the SPI in relation to the SDI on (A) 1 & (B) 3 month 
time scale for the Flint River Watershed (1970 to 2000). 
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Figure 65. Least-squares linear regression of the SPI in relation to the SDI on (A) 6 & (B) 9 month 
time scale for the Flint River Watershed (1970 to 2000). 
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Figure 66. Least-squares linear regression of the SPI in relation to the SDI on (A) 12 & (B) 
month time scale for the Flint River Watershed (1970 to 2000). 
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Figure 67. Least-squares linear regression of the SPI in relation to the SDI on (A) 1 & (B) 3 month 
time scale for the Middle Oconee River Watershed (1970 to 2000). 
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Figure 68. Least-squares linear regression of the SPI in relation to the SDI on (A) 6 & (B) 9 month 
time scale for the Middle Oconee River Watershed (1970 to 2000). 
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Figure 69. Least-squares linear regression of the SPI in relation to the SDI on (A) 12 & (B) 24 
month time scale for the Middle Oconee River Watershed (1970 to 2000). 
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Appendix D Metadata Sources 


Table 11. Summary of weather, streamflow, and ArcGIS Metadata. 


Description 

Data Name 

Data Type 

Topography 

n34w083-86, n35w083-86, & 
n34w083-86 

1 m DEM Raster 

Streamflow 

Surface-Water Daily Data 

Excel 

Weather 

Daily Summaries 

Excel 

2001 Land Use 

nlcd200 l_landcover 

Raster 

Streams 

NHD Best Resolution 

Shapefile - Line 

Chattooga Dams 

DamsChattooga 

Shapefile - Polygon 

Ecoregion Level III 

NA_CEC_Eco_Level3 

Shapefile - Polygon 

Ecoregion Level IV 

EcoLe veil V_U S 

Shapefile - Polygon 

Geology 

kbge 

Shapefile - Polygon 

Georgia Dams 

NID 

Shapefile - Polygon 

States 

State 

Shapefile - Polygon 

Topography 
Quadrant Names 

ned_larcsec_g 

Shapefile - Polygon 

Watersheds 

OW/WBDNAD83 

Shapefile - Polygon 


Data Location 

http://viewer.nationalmap.gov/viewer/ 

http ://waterdata. usgs. go v/n wis/s w 
https://www.ncdc.noaa.gov/data-access/land-based-station-data 
http ://www. mrlc. go v/nlcdO 1 data, php 
http://viewer.nationalmap.gov/viewer/ 

Created by author from aerial analysis. 
http://www.epa.gov/wed/pages/ecoregions/level_iii_iv.htm 
http://www.epa.gov/wed/pages/ecoregions/level_iii_iv.htm 
http://mrdata.usgs.gov/geology/us/ 
http://nid.usace.army.mil 
http://viewer.nationalmap.gov/viewer/ 
http://viewer.nationalmap.gov/viewer/ 

http://water.epa.gov/scitech/datait/tools/waters/services/mapping_services.cfm 















































